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We provide a systematic quantitative description of the structure of edge states and magnetosub- 
band evolution in hard wall quantum wires in the integer quantum Hall regime. Our calculations 
are based on the self-consistent Green's function technique where the electron- and spin interac- 
tions are included within the density functional theory in the local spin density approximation. We 
analyze the evolution of the magnetosubband structure as magnetic field varies and show that it 
exhibits different features as compared to the case of a smooth confinement. In particularly, in the 
hard- wall wire a deep and narrow triangular potential well (of the width of magnetic length Ib) is 
formed in the vicinity of the wire boundary. The wave functions are strongly localized in this well 
which leads to the increase of the electron density near the edges. Because of the presence of this 
well, the subbands start to depopulate from the central region of the wire and remain pinned in the 
well region until they are eventually pushed up by increasing magnetic field. We also demonstrate 
that the spin polarization of electron density as a function of magnetic field shows a pronounced 
double-loop pattern that can be related to the successive depopulation of the magnetosubbands. In 
contrast to the case of a smooth confinement, in hard- wall wires the compressible strips do not form 
in the vicinity of wire boundaries and spatial spin separation between spin-up and spin-down states 
near edges is absent. 

PACS numbers: 73.21.Hb, 73.43.-f, 73.23.Ad 



I. INTRODUCTION 

Recent advances in fabrication of low-dimensional 
structures allow one to create quantum wires with a hard- 
wall potential confinement. The available technologies 
include implantation-enhanced intcrdiffusion technique 1 
developed more than 20 years ago. Using this technique 
Prins et al& demonstrated a potential jump at a hetero- 
interface GaAs-AlGaAs over only 8 nm distance. The 
molecular beam epitaxy double-growth technique^ (of- 
ten referred to as a cleaved-edge overgrowth) since early 
1990-th has become one of the most widely- used tech- 
niques for fabrication of quantum wires^u^ and two- 
dimensional electron gases (2DEGs)£ with an essentially 
hard wall confinement with the atomic precision. Quan- 
tum wires with a steep confinement can also be fabricated 
by overgrowth on patterned GaAs(OOl) substrates using 
molecular beam epitaxy 8 -. 

For theoretical description of the quantum Hall effect 
in quantum wires, a concept of edge states is widely 
used^. In a naive one-electron picture a position of the 
edge states are determined by the intersection of the Lan- 
dau levels (bent by the bare potential) with the Fermi 
energy, and their width is given by a spatial extension 
of the wave function, which is of the order of the mag- 
netic length Ib = \J -^g- For a smooth electrostatic con- 
finement that varies monotonically throughout the cross- 
section of a wire, Chklovskii at ulri— have shown that elec- 
trostatic screening in strongly modifies the structure of 
the edge states giving rise to interchanging compressible 
and incompressible strips. The electrons populating the 
compressible strips screen the electric field, which leads 
to a metallic behavior when the electron density is re- 



distributed (compressed) to keep the potential constant. 
The neighboring compressible strips are separated from 
each other by insulator-like incompressible strips corre- 
sponding to the fully filed Landau levels with a constant 
electron density. 

A number of studies of quantum wires with a 
smooth confinement have been reported during the recent 
decade 11 ! 12 ! 13 . 14 : 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 addressing the prob- 
lem of electron-electron interaction beyond Chklovskii 
at al'sML electrostatic treatment. A particular atten- 
tion has been paid to spin polarization effects in the 
edge statea 11 ! 13 ! 16 ! 17 ! 19 ' 23 ! 24 . It has been demonstrated 
that the exchange and correlation interactions dramat- 
ically affect the edge state structure in quantum wires 
bringing about qualitatively new features in comparison 
to a widely used model of spinless electrons. These in- 
clude spatial spin polarization of the edge statesi 3 * 2 ^, pro- 
nounced l/£?-pcriodic spin polarization of the electron 
density^ 3 ., modification and even suppression of the com- 
pressible stripsSi and others. It should be stressed that 
all the above-mentioned studies addressed the case of a 
soft confinement corresponding to e.g. a gate-induced 
depletion when the Borh radius is much smaller that the 
depletion length. In fact, Huber et al. have recently pre- 
sented experimental evidence that widely used concept 
of compressible/incompressible strips 1 ^ does not apply 
to the case of a sharp-edge 2DEG. At the same time 
the rigorous theory for edge-state structure in hard-wall 
quantum wires accounting for electron-electron interac- 
tion and spin effects has not been reported yet. Such a 
theory is obviously required for a detailed analysis of re- 
cent experiments on cleaved-edge overgrown sharp-edge 
wires and 2DEGo 2 i 3 i 4 i 5 i 6 i 7 i 8 . 

Motivated by the above-mentioned experimental stud- 
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ics, in this paper wc present a detailed theory of magneto- 
subband and edge state structure in quantum wires with 
a hard wall confinement taking into account electron- 
electron interaction including exchange and interaction 
effects. We employ an efficient numerical tool based on 
the Green's function technique for self-consistent solu- 
tion of the Schrodinger equation in the framework of the 
density functional theory (DFT) in the local spin density 
approximation (LSDA)S. The choice of DFT+LSDA 
for description of many-electron effects is motivated, 
on one hand, by its efficiency in practical implementa- 
tion within a standard Kohn-Sham formalism^, and, on 
the other hand, by an excellent agreement between the 
DFT+LSDA and the exact diagonalizationSi and the 
variational Monte-Carlo calculations^*^ performed for 
few-electron quantum dots. We will demonstrate below 
that edge state structure of the hard wall quantum wire is 
qualitatively different from that of the soft-wall wire. We 
will discuss how the spin-resolved subband structure, the 
current densities, the confining potentials, as well as the 
spin polarization in the hard wall quantum wire evolve 
when an applied magnetic field varies. 

The paper is organized as follow. In Sec. II we 
present a formulation of the problem, where wc define 
the geometry of the system at hand and outline the self- 
consistent Kohn-Sham scheme within the DFT+LSDA 
approximation. In Sec. Ill we present our results for a 
hard wall quantum wire calculated within Hartree and 
DFT+LSDA approximations, where we distinguish cases 
of wide and narrow wires. Section IV contains our con- 
clusions. 
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FIG. 1: (Color online), (a) A schematic illustration of a 
cleaved-edge overgrown quantum wire and (b) a correspond- 
ing hard-wall confinement potential. 



framework of the Kohn-Sham density functional theory 
reads2&, 



V: ff (y)=V H (y) + V: x (y), 



(3) 



where Vn(y) is the Hartree potential due to the elec- 
tron density n(y) = ^ a n a (y) (including the mirror 
charges)^, 



II. MODEL 

We consider a quantum wire which is infinitely long in 
the x-dircction and is confined by a hard-wall potential 
in the y-direction, see Fig. ^ 

The magnetic field is applied perpendicular to the xy- 
plane. We set the Fermi energy Ep = 0. A bottom of the 
confining potential is flat and situated at E = Vq. We 
limit ourself to a typical case when only one subband is 
occupied in the transverse z-directioni such that electron 
motion is confined to the xy-planc. The Hamiltonian of 
the wire reads H = Y2cr H c r , 

H° = H + V + Vf ff (y) + g^Ba, (1) 

where Hq is the kinetic energy in the Landau gauge, 

h 2 



Hn = — 



2m* 




where a = ±^ describes spin- up and spin-down states, 1, 
I, and to* = 0.067m e is the GaAs effective mass. The last 
term in Eq. (1) accounts for Zeeman energy where [ib = 

tt~~ — is the Bohr magneton, and the bulk q factor of GaAs 

2m e ° ' & 

is g = —0.44. The effective potential, V e ff(y) within the 
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with 2b being the distance from the electron gas to the 
mirror charges (we choose b=60 nm). For the exchange 
and correlation potential V xc {y) we utilize the widely 
used parameterization of Tanatar and Cerperly2£ (see 
Ref. |22| for explicit expressions for V xc (y))- This pa- 
rameterization is valid for magnetic fields corresponding 
to the filling factor v > 1, which sets the limit of appli- 
cability of our results. The spin-resolved electron density 
reads 



n a (y) 



^ J dEG a (y,y,E)f FD (E-E F ), 



(5) 



where G cr (y, y, E) is the retarded Green's function corre- 
sponding to the Hamiltonian and fFD{E — Ep) is the 
Fermi-Dirac distribution function. The Green's function 
of the wire, the electron and current densities are cal- 
culated self-consistently using the technique described in 
detail in Ref. 

The current density for a mode a is calculated as 23 
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with v° a and j„(y, E) being respectively the group veloc- 
ity and the quantum-mechanical particle current density 
for the state a at the energy E, and V being the applied 
voltage. 

We also calculate a thermodynamical density of states 
(TDOS) defined according to2±^ 

TDOS* = J dEp°{E) {- dfFDi g E EP) ) , (7) 

where the spin-resolved density of states p a {E) is given 
by the Green function^, 

p°(E) = -~%[dyG°(y,y,E). (8) 

The TDOS reflects a structure of the magnetosubbands 
near the Fermi energy and it can be accessible via 
magneto-capacitance^ or magnetoresistance^ measure- 
ments. Indeed, a compressible strip corresponds to a 
flat (dispersionless) subband pinned at Ep. In this case 
p a {E) is high at E w Ep and such the subband strongly 
contributes to TDOS. In contrast, in an incompress- 
ible strip, subbands are far away from Ef and do not 
contribute to TDOS. Thus TDOS is proportional to 
the area of the compressible strips. This area is maximal 
when the strip if formed in the middle of a quantum wire. 
In this case the backscattering between opposite propa- 
gating states is maximal, which corresponds to peaks in 
the longitudinal resistance R xx (seen as the Shubnikov- 
De Haas oscillations)2£i2£i2I. In magneto-capacitance 
experiments^^ the compressible strips are viewed as 
capacitor plates and therefore the measured magnetoca- 
pacitance is related to the width of these strips. Thus 
the peaks in the TDOS are manifest themselves in both 
R xx and capacitance peaks. 

III. RESULTS AND DISCUSSION 

In what follows we shall distinguish between cases of 
a wide quantum wire whose half-width ^ exceeds the 
magnetic length / b , and a narrow wire with a width ^ < 
l B . 

A. Wide hard wall quantum wire ^ > Ib 

Let us consider a hard wall quantum wire of the width 
w = 300 and Vq = —0.1 eV. With these parameters the 
wire has N ~ 20 spin-resolved occupied subbands at zero 
magnetic field, and the sheet electron density in its center 
is ri2D ~ 1-5 • 10 15 m~ 2 (as calculated self-consistently in 
both Hartree and DFT approximations). 

a. Hartree approximation We start our analysis of 
the edge state- and magnetosubband structure from the 
case of the Hartree approximation (when the exchange 
and correlation interactions are not included in the ef- 
fective potential). The Hartree approximation gives 



the structure of the compressible/incompressible strips 
which serves as a basis for understanding of the ef- 
fect of the exchange and correlation within the DFT 
approximation2ii24. 

Figure Ola) shows the ID electron density n a D = 
J n a (y)dy for the spin- up and spin-down electrons in the 
quantum wire. The pronounced feature of this depen- 
dence is a characteristic loop pattern of the charge den- 

x _ l 

sity polarization, P n = n ) D n \ D , see Fig. EJb) . Fig- 

ure |2l also indicates a number of magnetosubbands N 
populated at a given B. The number of subbands is 
always even such that spin-up and spin-down subbands 
depopulate practically simultaneously. This is because 
the spin polarization within the Hartree approximation 
is driven by Zeeman splitting only, which is small in the 
field interval under consideration. A comparison of Figs. 
|2{a),(c),(e) demonstrates that the spin polarization as 
well as the TDOS are directly related to the magneto- 
subband structure. Note that a similar loop-like behavior 
of the spin polarization is also characteristic for a split- 
gate wire with a smooth confinement^. For the latter 
case the polarization calculated in the Hartree approxi- 
mation drops practically to zero when the subbands de- 
populate (see Fig. 4 in Ref. I23T) . In contrast, in the 
case of the hard wall confinement, the polarization loops 
exhibit more complicated pattern: the polarization does 
not drop to zero when the subbands depopulate, and, in 
addition, the polarization curves show a double loop-like 
pattern with an additional minimum (e.g. at B ss 1.5T, 
3T in Fig. [U(a),(c)). In order to understand the origin 
of this behaviour let us analyze the evolution of the sub- 
band structure as the applied magnetic field varies. Let 
us concentrate at the field interval 1.65 T < B < 3.5 T 
when the subband number N = 4. 

Figure E2b) shows the spatially resolved difference in 
the electron density rv(y) — n^(y) as a function of B. 
The electron density is mostly polarized in the inner re- 
gion of the quantum wire. For certain ranges of magnetic 
fields the electron density shows a strong polarization in 
the boundary regions, which are separated from the po- 
larized inner region by wide unpolarizcd strips (e.g. for 
3T < B < 3.5T). We will show below that this feature 
reflects the peculiarities of the magnetosubband struc- 
ture for the case of the hard wall confinement. Figure 
I3c) shows the electron density profiles (local filling fac- 
tors) v{y) = n(y)/riB {ns = eB/h), the current den- 
sities J a (y) and the magnetosubband structure for the 
magnetic field B = 1.8 T. At this field a wide com- 
pressible strip due to electrons belonging to the subbands 
N = 3, 4 is formed in the middle of the wire. (Following 
Suzuki and Ando^S we define the width of the compress- 
ible strips within the energy window \E — Ep\ < lixkT 
corresponding to the partial occupation of the subbands 
when fpr) < 1; this energy window is indicated in Fig. 
(c)). Partial subband occupation combined with Zeeman 
splitting of energy levels results in different population 
for spin-up and spin-down electrons (i.e. in the spin po- 
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FIG. 2: (Color online). (a),(b) One-dimensional charge density for the spin-up and spin-down electrons, n\ D , n[ D ; (c),(d) the 



spin polarization of the charge density, P n 



T _ I 
1 1,D "lD 



(g),(h) the TDOS for spin- up and spin-down electrons and the total 



TDOS within the Hartree approximation and the DFT approximation (first and second columns, respectively). The number 
of subband is indicated in (a),(b). Arrows in (c) and (d) indicate the magnetic field corresponding to the magnetosubband 
structure shown in Figs. [^]and[4] The width of the wire is w — 300 nm and the depth is Vo = —0.1 eV. Temperature T = 1 K. 



larization of the electron density). 

Close to the wire edges the total potential exhibits a 
narrow and deep triangular well. The formation of the 
triangular well is also reflected in the structure of the 
magnetosubbands that show triangular wells near the 
wire edges. Presence of these triangular wells is a dis- 
tinctive feature of the hard- wall confinement (it is absent 
for the case of a smooth confinement in the split-gate 
wiresi2i2L2£2i) . The wave functions for all subbands are 
strongly localized in these wells, with the extension of 
the wave functions being of the order of the magnetic 
length Ib- Because of steepness of the potential walls, the 
wave functions are not able to screen the confining poten- 
tial, and compressible strips can not form near the wire 
boundary. This is in a stark contrast to the case of a split- 
gate wire where the compressible strips near edges are 
formed for a sufficiently smooth confincmentpiS*22i2&2i. 
The electron density near the wire boundaries does not 
show any spin polarization. This is because the bottom 
of the potential well lies far below the Fermi energy. As 
a result, both spin-up and spin-down states localized in 
the quantum well are completely filled (/fd = 1) and 
the spin polarization is absent. 

When a magnetic field increases the compressible strip 



in the middle of the wire widens. This is accompanied 
by increase of both the spin polarization and the TDOS 
as shown in Figs. Etc),(e). At B = 2.3T the polarization 
reaches maximum P n = 3 % which corresponds to the 
maximum width of the compressible strip in the central 
part of the wire, see Fig. IHId). With further increase 
of the magnetic field 3rd and 4th subbands in the cen- 
tral part of the wire are pushed up, see Fig. Efe). Their 
population decreases according to the Fermi-Dirac dis- 
tribution and, consequently, the spin polarization dimin- 
ishes. At the same time, fully occupied parts of 3rd and 
4th subbands (forming a triangular well near the wire 
boundaries) arc pushed up and got pinned at the Fermi 
energy. This is accompanied by a formation of a poten- 
tial barrier at the distance of the wave function extent 
~ Ib from the wire edges, see Fig. E^e). The whole area 
occupied by subbands 3 and 4 becomes divided by non- 
populated region within the barrier where the subbands 
lie above Ep (i. e. fpp, = 0). 

When a magnetic field slightly increases from B = 
2.8T to B = 3.0T the magnetosubband structure under- 
goes significant changes. A middle part of the 3rd and 
4th subbands is abruptly pushed up in energy. The in- 
compressible strip emerges here due to 1st and 2nd fully 



5 



0- = 



(a) 


6 










r 



(10 m ) 



10 




1 


2 3 


4 




B(T) 







>14 




100 
y (nra) 



100 
y (nm) 



100 
y (nm) 



100 

y (nm) 



FIG. 3: (Color online), (a) Spin polarization of the charge density as a function of B calculated within the Hartree approximation 
(the same as Fig. EJc)). (b) Spatially resolved difference in the electron density n' (y) — n^(y). (c)-(g) The subband structure 
for magnetic fields indicated by arrows in (a). Upper panel: electron density profiles (local filling factors) v(y) = n(y)/riB 
for spin-up and spin-down electrons; middle panel: the current density distribution for spin-up and spin-down electrons; lower 
panel: magnetosubband structure for spin-up and spin-down electrons. Fat solid and dashed lines indicate the total confining 
potential for respectively spin-up and spin-down electrons. The width of the wire is a — 300 nm and depth is Vq = —0.1 eV. 
Temperature T=l K. 



occupied subband lying well below Ep, Fig. EHf). As a 
result the spin polarization decreases and the fist polar- 
ization loop closes down at B w 3T, see^a). Note that 
P n does not drop to zero because of a finite polariza- 
tion at the boundaries where the 3rd and 4th subband 
bottoms arc still pinned at the Fermi energy, see Fig. 
EIb),(f). As magnetic field increases the second polariza- 
tion loop starts to form at B « 3T due to 1st and 2nd 
subbands that get pinned to Ef in the middle of the wire 
(Fig. EUg)). I n addition, 3rd and 4th subbands that are 
pinned to Ef near the wire boundaries also contribute 
to spin polarization. These subbands become completely 
depopulated at B = 3.5 T. Further increase of the mag- 
netic field causes the compressible strip in the middle to 
widen. The spin polarization P n grows linearly until the 
second subband becomes depopulated. 

Note that the above scenario of the subband depopu- 
lation in quantum wires with a hard wall confinement is 
qualitatively different from that one of the smooth con- 
finement. In the former case, because of the presence of 
the deep triangular well near the wire boundaries, the 
subbands start to depopulate from the central region of 
the wire and remain pinned in the well region until they 
are eventually pushed up by magnetic field. In contrast, 
in the case of a smooth confinement, the subband always 
depopulate from the edges, such that a compressible strip 
in the middle of the wire gradually decreases until it com- 
pletely disappears when the whole subband is pushed up 
above the Fermi energy22i24. 



The spatial current distribution stays practically the 
same throughout the magnetosubband evolution, see the 
central panel in Figs. Etc)-(g). This is due to a strong 
localization of electrons in the triangular potential well. 
The spatial spin separation between spin-up and spin- 
down states is always equal to zero, which is also the case 
for a split-gate wire in the Hartree approximation^^. 

Finally, within the Hartree approximation the TDOS 
shows a behavior similar to the spin polarization of the 
electron density P n , compare Fig. |5fe) andHJc). This 
is because the spin polarization is primarily caused by 
electrons in the compressible strips, and the TDOS, as 
discussed in the previous section, is proportional to the 
width of these strips. 

b. DFT approximation The exchange and correla- 
tion interactions bring qualitatively new features to the 
magnetosubband structure in comparison to the Hartree 
approximation. Figures|2fb),(d),(f) show the ID electron 
density, the number of subbands, the spin polarization 
and the TDOS calculated within DFT approximation. 
There are several major differences in comparison to the 
Hartree case. First, the spin polarization of the electron 
density also shows a pronounced loop pattern. However, 
for a given magnetic field the spin polarization in the 
quantum wire calculated on the basis of the DFT ap- 
proximation is much higher in comparison to the Hartree 
approximation (by a factor 5-10). Second, the exchange 
interaction lifts subband degeneration, such that the sub- 
bands depopulate one by one. Third, the TDOS reveals 
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FIG. 4: (Color online), (a) Spin polarization of the charge density as a function of B calculated within the DFT approximation 
(similar to Fig. EJd)). (b) Spatially resolved difference in the electron density n T (y) — n^(y). (c)-(g) The subband structure 
for magnetic fields indicated by arrows in (a). Upper panel: electron density profiles (local filling factors) v(y) = n{y)/riB 
for spin-up and spin-down electrons; middle panel: the current density distribution for spin-up and spin-down electrons; lower 
panel: magnetosubband structure for spin-up and spin-down electrons. Fat solid and dashed lines indicate the total confining 
potential for respectively spin-up and spin-down electrons. The width of the wire is a = 300 nm and depth is Vq = —0.1 eV. 
Temperature T — IK. 



peaks which are attributed to different spin species. 

Before we proceed to analysis of the magnetosubband 
structure within the DFT approximation, it is instru- 
mental to outline the effect of the exchange interaction 
on the subband spin splitting. Within the Hartree ap- 
proximation the subbands are practicably degenerate be- 
cause the Zeeman splitting is very small in the magnetic 
field interval under investigation. In contrast, the ex- 
change interaction included within the DFT approxima- 
tion causes the separation of the subbands which mag- 
nitude can be comparable to the Landau level spacing 
Huj. Indeed, the exchange potential for spin-up electrons 
depends on the density of spin-down electrons and vice 
verga 23 i 25 4 30^ j n compressible region the subbands 
are only partially filled (because jpD < 1 in the the win- 
dow \E — Ep\ < 2nkT) 1 and, therefore, the population 
of the spin-up and spin-down subbands can be differ- 
ent. In the DFT calculation, this population difference 
(triggered by Zeeman splitting) is strongly enhanced by 
the exchange interaction leading to different effective po- 
tentials for spin-up and spin-down electrons and even- 
tually to the subband spin splitting. Below the Fermi 
energy E < Ep — 2nkT the subbands remain degenerate 
because they are fully occupied (/fd = 1)- As a re- 
sult, the corresponding spin-up and spin-down densities 
are the same, hence the exchange and correlation poten- 
tials for the spin-up and spin-down electrons are equal, 
Vl(y) = V x {(y). 

In order to understand the effect of the exchange- 



correlation interactions on evolution of the magnetosub- 
band structure, let us concentrate on the same field in- 
terval as discussed in the case of the Hartree approxima- 
tion, 1.8 T < B < 3.7 T. A comparison between Fig. 0] 
and Fig. [3] demonstrates that evolution of the magne- 
tosubband structure calculated within the DFT approx- 
imation follows the same general pattern as for the case 
of the Hartree approximation. In particularly, a deep 
triangular well near the wire boundary develops in the 
total confining potential for both spin-up and spin-down 
electrons. The wave functions are strongly localized in 
this well. As a result, similarly to the Hartree case, the 
depopulation of the subbands starts from the central re- 
gion of the wire. The subbands remain pinned in the 
well region until they are eventually pushed up by mag- 
netic field. The major difference from the Hartree case 
is that Hartree subbands are practically degenerated and 
depopulate together, whereas this degeneracy is lifted by 
the exchange interaction such that DFT subbands de- 
populate one by one. Indeed, Figs. 0] (c),(d) showing 
consecutive depopulation of the subbands 4 and 3 in the 
central region of the wire can be compared with the cor- 
responding evolution of the Hartree subbands in Figs. 
(c),(d). When the magnetic field increases further, 3rd 
subband bends upward in the vicinity of the triangular 
well, compare Fig. 0](e) and Fig. |3(e). When magnetic 
field reaches B w 2.7T, 4th spin-down subband becomes 
completely depopulated and 3rd spin-up subband is oc- 
cupied mostly in the region of the triangular well near 
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the wire boundary, see Fig. 0(f). This leads to a strong 
spin polarization near the boundary which is manifest it- 
self in the additional loop of the polarization (see Fig. [3 
(b), 2.7T< B <3.2T). Note that this loop is absent in 
the Hartree calculations because both 3rd and 4th sub- 
bands are occupied in the well region, such that the spin 
splitting between them is small (see Fig. 01(f))- Finally, 
3rd subband becomes fully depopulated in the central re- 
gion, and a compressible strip starts to form there due 
to 2nd subband that is pushed upwards, compare Figs. 
11(g) and Fig. (g). 

Note that similarly to the case of the Hartree approx- 
imation, the evolution of the magnctosubband structure 
within the DFT approximation described above qualita- 
tively holds for all other polarization loops. 

We also stress that in contrast to the case of a smooth 
confinemenljSiiSi, in hard-wall wires the compressible 
strips do not form in the vicinity of wire boundaries and 
a spatial spin separation between spin-up and spin-down 
states near edges is absent. 

Oscillations of the TDOS calculated within the DFT 
approximation shows that neighboring peaks belong to 
different spin species (Fig. Eff)). In contrast, the Hartree 
approximation shows that each single peak includes equal 
contributions from both species (Fig. Gle)). It is in- 
teresting to note that the oscillation of the TDOS do 
not exactly correspond to the subband depopulation. In- 
stead, they reflect formation of the compressible strip in 
the middle of the wire due to spin-up and spin-down elec- 
trons which is not directly related to the subband depop- 
ulation (which takes place in the region of the triangular 
well near the wire edge). 

To conclude this section we note that we analyzed the 
magnctosubband structure for a representative sharp- 
edge quantum wire of 300 nm width. It is important 
to stress that all the conclusions presented above (i.e. 
the scenario of magnctosubband depopulation and the 
structure of the edge states near the wire boundary) hold 
for an arbitrary sharp-edged quantum wire provided its 
length is sufficiently larger than the magnetic length Ib- 
In particulary, our results can be applied to analysis of an 
cpitaxially overgrown cleaved edge semi-infinite structure 
similar to that one studied in Ref. 0- 



B. Narrow hard wall quantum wire ^ < Ib 

Let us now concentrate on the case of a narrow wire 
whose half-width is comparable to the magnetic length. 
For our analysis we choose the wire of the width w=50 nm 
and Vo = —0.2 eV. With these parameters the electron 
density at the center of the wire is n2D ~ 6 • 10 15 m -2 
and the number of spin-resolved subbands is N = 6 for 
B = T. 

Figures EJa) and (b) show respectively the ID charge 
density and the polarization for spin-up and spin-down 
electrons calculated within the DFT approximation. Let 
us concentrate on the field interval 7 < B < 12, when 



a number of subbands 3 < N < 4. In this interval the 
spin polarization shows a pronounced single-loop pattern. 
This is in contrast to the case of a wide wire that exhibits 
a double-loop pattern (see Figs. |21 (a),(b)), where the 
first loop corresponds to the subband depopulation in the 
middle of the wire, whereas the second loops corresponds 
to the subband depopulation in the deep triangular well 
near the boundary. Note that the width of the this well is 
of the order of the extension of the wave function given 
by the magnetic length l B . This explain a single- loop 
structure of the polarization curve for the case of a nar- 
row wire ^ < Ib- Indeed, in this case the extension of 
the triangular well is comparable to the half- width of the 
wire, such that the well extends in the middle region and 
there is no separate depopulation for the inner and outer 
regions of the wire. 

The above features of the narrow wire can be clearly 
traced in the evolution of the magnetosubbands, see Fig. 
El When 6.5 T < B < 8.5 T 3rd and 4th subbands in 
the middle of the wire are located beneath Ep — 2irkT 
and are thus fully occupied. This corresponds to the for- 
mation of the incompressible strip in the middle of the 
wire such that the charge densities of spin-up and spin- 
down electrons are equal (i.e. the spin polarization is 
zero). At B — 8.5T 4th subband reaches E F - 2irkT 
and thus becomes partially occupied. As a result, the 
exchange interactions generates spin splitting, and the 
compressible strip due to spin-down electrons belonging 
to 4th subband starts to form in the middle of the wire. 
Spin polarization grows rapidly until it reaches its max- 
imum P n = 22 %. At this moment 4th subband de- 
populates and the corresponding compressible strip dis- 
appears. When magnetic field is increased only slightly, 
3rd subband is raised to Ep — 2irkT and the compressible 
strip due to spin-up electrons forms in the middle of the 
wire. Note that formation and disappearance of the com- 
pressible strips due to spin-up and spin-down electrons is 
clearly reflected in the TDOS, see Fig. [51(c) which shows 
peaks belonging to different spin species. With further 
increase of B the spin polarization decreases linearly until 
it vanishes when 3rd subband fully depopulates. 

Magnctosubband evolution calculated within the 
Hartree approximation (not shown) qualitatively resem- 
bles evolution for the DFT case. In particular, the spin 
density polarization follows the same behavior reaching 
the maximum value P n = 10 in the interval 3 < N < 4. 
The similarity between the Hartree and DFT approxima- 
tions is because of a large Zeeman term for magnetic field 
intervals under consideration which causes a relatively 
strong Zeeman splitting in the Hartree approximation. 



IV. CONCLUSION 

We provide a systematic quantitative description of the 
structure of the edge states and magnetosubband evolu- 
tion in hard wall quantum wires in the integer quan- 
tum Hall regime. Our calculations are based on the 



8 





10 20 10 20 10 20 
y (nm) y (nm) y (nm) 

FIG. 5: (Color online). ID charge density for spin-up and spin-down electrons (a), the charge spin polarization (b), the TDOS 
for spin-up and spin-down electrons, total TDOS (c) as a function of B calculated within the DFT approximation for a narrow 
wire, (d)-(f) The subband structure for magnetic fields indicated in (b). Upper panel: electron density profiles (local filling 
factors) v(y) — n(j/)/ns; middle panel: the current density distribution; lower panel: magnetosubband structure for spin-up 
and spin-down electrons. Fat solid and dashed lines indicate the total confining potential for respectively spin-up and spin-down 
electrons. The width of the wire is a = 50 nm and depth is Vo = —0.2 eV. Temperature T = IK. 



self-consistent Green's function technique^ where the 
electron- and spin interactions are included within the 
density functional theory in the local spin density ap- 
proximation. Our main findings can be summarized as 
follows. 

1) The magnetosubband structure and the density dis- 
tribution in the hard-wall quantum wire is qualitatively 
different from that one with a smooth electrostatic con- 
finement. In particularly, in the hard- wall wire a deep 
triangular potential well of the width ~ Ib is formed in 
the vicinity of the wire boundary. The wave functions are 
strongly localized in this well which leads to the increase 
of the electron density near the edges. 

2) Because of the presence of the deep triangular well 
near the wire boundaries, the subbands start to depopu- 
late from the central region of the wire and remain pinned 
in the well region until they are eventually pushed up by 
an increasing magnetic field. This is in contrast to the 
case of a smooth confinement where depopulation of the 



subbands starts from the edges and extends towards the 
wire center as the magnetic field increases. 

3) The spin polarization of electron density as a func- 
tion of magnetic field shows a pronounced double-loop 
pattern that can be related to the successive depopula- 
tion of the magnctosubbands. 

4) In contrast to the case of a smooth confinement, in 
the hard-wall wires the compressible strips do not form 
in the vicinity of wire boundaries and a spatial spin sep- 
aration between spin-up and spin-down states near the 
edges is absent. 
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